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ABSTRACT 



We use dimensional regularization to calculate the Ois^^ expansion of all scalar one- loop one-, 
two-, three- and four-point integrals that are needed in the calculation of hadronic heavy quark 
production. The Laurent series up to 0{£^^ is needed as input to that part of the NNLO cor- 
rections to heavy flavor production at hadron colliders where the one-loop integrals appear in the 
loop-by-loop contributions. The four-point integrals are the most complicated. The 0(s'^^ ex- 
pansion of the three- and four-point integrals contains in general polylogarithms up to Li4 and 
functions related to multiple polylogarithms of maximal weight and depth four. 
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I. INTRODUCTION 



The full next-to-leading order (NLO) corrections to the hadroproduction of heavy fla- 
vors have been completed in the late eighties [1,2]. They have raised the leading order (LO) 
estimates [3] but several initial analysis' showed a serious disagreement with experimental 
results [4, 5]. Recently the situation has considerably improved in that a more refined 
NLO analysis (due to considerably more precise experimental input for the 6-quark frag- 
mentation function as well as other QCD parameters) now shows signs of rapprochement 
between theory and the new experimental data (see [6] and references therein for the new 
CDF measurements). However, the NLO predictions are still slightly below the experi- 
mental numbers. Moreover, the theoretical NLO predictions suffer from the usual large 
uncertainty resulting from the freedom in the choice of renormalization and factorization 
scales of perturbative QCD. In this hght there are hopes that a next-to- next-to-leading 
order (NNLO) calculation will bring theoretical predictions even closer to the experimental 
data. Also, the dependence on the factorization and renormalization scales of the physical 
process is expected to be greatly reduced at NNLO. This would reduce the theoretical un- 
certainty and therefore make the comparison between theory and experiment much more 
significant. 

In Fig. 1 we show one generic diagram each for the four classes of gluon-induced contri- 
butions that need to be calculated for the NNLO corrections to hadroproduction of heavy 
flavors. They involve the two-loop contribution (la), the loop-by-loop contribution (lb), 
the one-loop gluon emission contribution (Ic) and, finally, the two gluon emission contri- 
bution (Id). A similar classification holds for the quark-induced contributions. 




c) d) 



FIG. 1. Exemplary gluon fusion diagrams for the NNLO calculation of heavy hadron production. 

In this paper we concentrate on the loop-by-loop contribution Fig. lb. Specifically, 
working in the framework of the dimensional regularization scheme [7], we shall present 
0{e'^) results on all scalar one- loop one-, two-, three- and four-point integrals that are 
needed in the calculation of hadronic heavy fiavour production. We generate the coefficients 
of the e-expansion in a rather direct way. Let us briefiy describe our procedure. We 
introduce Feynman parameter representations for each of the two-, three- and four-point 



1 



integrals. The two-point case is straightforward. In case of the three- and four-point 

integrals we generally integrate over one and two Fcynman parameters, respectively, keeping 
the full ^-dependence of the result. Before doing the last Feynman parameter integration 
we expand the respective integrands up to 0{e'^) and then integrate the expanded integrand 
term by term. 

Because the one-loop integrals exhibit ultraviolet (UV) and infrared (IR)/coUinear 
(or mass (M)) singularities up to 0{e~'^) one needs to know the one-loop integrals up 
to because the one-loop contributions appear in product form in the loop-by-loop 

contributions^. It is clear that the spin algebra and the Passarino-Veltman decomposi- 
tion of tensor integrals in the one-loop contributions also have to be done up to 
This task will be left to a companion paper where we present complete results on the one- 
loop amplitudes up to 0{e'^) including spin algebra and Passarino-Veltman decomposition 
effects. 

The general case of massive one-loop integrals was studied some time ago [9] , where a 
general one-loop A^-point integral was expressed in terms of hypergeometric functions of 
several variables. Recently, there have been a number of papers where the authors took 
a more general attitude to calculate the e-expansion of massive one-loop integrals. They 
write down general representations of the e-expansion of one-loop integrals for general 
kinematic configurations. We have attempted to compare our results with the results of 
the more general approaches whenever possible. In papers [10, 11] the all-order £-expansion 
of one-loop two-point and of certain three-point functions was done explicitly by expanding 
the relevant hypergeometric functions. One-fold integral representations for general three- 
and four-point functions, as well as ways to get expansion terms of order e for 3-point 
functions, were worked out in a recent paper [12]. Publications [9, 11, 12] also contain a 
comprehensive hst of references on the subject. 

However, in general, the required ^-expansion (including £^-terms) is not readily avail- 
able for all the integrals needed in the hadronic heavy flavour production process. Also, 
the analytic continuation of the above mentioned hypergeometric functions in [9, 12] to 
the appropriate kinematical regions of validity is not always possible. This mainly con- 
cerns the four-point functions. In addition, it is more convenient to present results for 
the ^-expansion in terms of simpler special functions, in the form convenient for numerical 
evaluation. And finally, collecting together all the necessary scalar integrals needed for the 
derivation of tensor integrals entering the loop-by-loop contribution constitutes a first step 
in the difficult task of obtaining the NNLO corrections to heavy flavor hadroproduction 
cross section. 

In our notation we shall remain very close to the notation introduced and used in [2,13]. 
In particular, we use dimensional regularization working in n = 4 — 2^ dimensions, as e.g. 
in [13]. For the calculation of the NNLO virtual corrections to hadroproduction of heavy 

"'^In a more general setting the Laurent-series expansion of the scalar integrals is needed if the integration- 
by-parts technique [8] is employed. The reason is that the solution of the recursion relations induced by 
the integration- by-parts technique can bring in negative powers of e. 
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flavors one needs the same set of scalar master integrals as given in the Appendix A of 
[2] (the relevant set of master integrals is listed in Table 1). However, as explained above, 
knowledge of their singular and finite terms is not sufficient for the calculation of NNLO 
loop-by-loop corrections. For that purpose one needs to know the one-loop integrals up to 
0{e^) including also their imaginary parts which equally well contribute to the modulus 
squared of the one-loop amplitudes. The imaginary parts of the one-loop integrals are 
really needed only up to 0{e) since the highest singularity of the imaginary parts is only 
0{e~^) compared to 0{e^'^) for the real parts. We have nevertheless decided to include 
0{e^) results also for the imaginary parts which may be of interest in other applications. 
Consequently, in this paper we present the relevant expressions for all scalar integrals 
needed in the calculation of the NNLO loop-by-loop corrections to hadroproduction of 
heavy ffavors. For reasons of comprehensiveness we have decided to include also the singular 
and finite (i.e. 0{e^)) parts of the scalar integrals in our presentation. They agree with 
the results of the real contributions presented in [2] . 

A comment on the length of the formula expressions in our paper is appropriate. The 
untreated computer output of the integrations is generally quite lengthy. The hard work is 
to simplify these expressions. We have written semi-automatic computer codes that achieve 
the simplifications using known identities among polylogarithms and using a number of 
identities for the L-functions introduced in this paper which are derived in an Appendix. 

The paper is organized as follows. In Sec. II we deal with one- and two-point functions. 
Sec. Ill contains our results on three-point functions, while in Sec. IV we present our re- 
sults for the ^-expansion of the four-point functions. In Sec. V we give our summary and 
conclusions. We collect some technical material in three Appendixes. In Appendix A we 
discuss the Taylor series expansion around = w? for the self-energy insertion two-point 
function which is needed for the calculation of the heavy quark wave function renormaliza- 
tion constant. In Appendix B we define multiple polylogarithms and demonstrate that our 
analytical results can all be expressed in terms of multiple polylogarithms. In Appendix C 
we derive a number of identities for the so-called L-functions introduced in the main text. 
A judicious use of these identities has allowed us to considerably reduce the length of our 
final analytical results for the three- and four-point functions. 



Working in n — A — 2s dimensions, the expansion for this one-point scalar integral can be 
written in the general form: 



II. ONE- AND TWO-POINT FUNCTIONS 



We start with the one-loop one-point function which is defined by 





A{m) = iCe{rn^) 



e{l-e) 



iC,{m^)m^ (i + l + £ + £2 + 0{e'') 



(2.2) 
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Table 1: List of one-, two-, three- and four-point massive one-loop functions calculated in 
this paper up to 0{e'^). 





Nomenclature of [2\ 


Our nomenclature 


Novelty 


Comments 


1-point 


A{m) 


A 


— 


Re 


2-point 


B{pa -P2,0,m) 


Bi 


— 


Re 




B{p3 + pA-,m,m) 


B2 


— 


Re, Im 




B{pA,0,m) 


B3 


— 


Re 




B{p2,m,m) 


B4 


— 


Re 




5(P3+P4,0,0) 


B5 




Re, Im 


3-point 


C(P4,P3,0,m, 0) 


^\ 


new 


Re, Im 




C(P4, -P2,0,m,m) 


C2 


new 


Re 




C{-p2,P4,0,0,m) 


C3 




Re 




<^(-P2, -Pi, 0,0,0) 


C4 




Re, Im 




C(-P2, -pi,m,m,m) 


C5 




Re, Im 




C{p3,p4,m,0,m) 


Ce 




Re, Im 


4-point 


D{p4, -p2, -pi, 0, m, m, m) 


Di 


new 


Re, Im 




^(-P2,P4,P3,0,0,m,0) 


D2 


new 


Re, Im 




D{-P2,P4, -pi,0,Q,m,m) 


Ds 


new 


Re 



where m is the internal loop mass and where we have defined 



r(l + £) Unfi^ 
(47r)2 y 



■ (2.3) 



The one-loop two-point functions are defined by [2] 

/c?^^ 1 
Tr^iT~2 2VT^ \2 2T' (2-4) 

where the [i = 1, 2) can be either m or 0. In the denominators of the relevant functions 
we always imply the "causal" +i6 prescription to deal with singularities in pseudo- Euclidean 
space. 

In what follows, we will always present our results for the scalar functions separately 
for the real and imaginary contributions. We introduce the Mandelstam-type variables 

s = {pi+p2f, t = T -m^ = {pi-p^Y -m^ , u = U -w? = {P2-P3f -w? , {2.b) 

with the kinematical condition on external momenta being pi+p2 = P3+P4 (i-e. s+t+u = 0) 
and the on-shell conditions are pf = p?2 = 0. p\ = p\ = 'rn?- Note that the variables t and 
u defined in (2.5) are not the usual Mandelstam variables. 

There are altogether five different two-point scalar functions B^ [i — 1, 2, 5) needed 
for hadronic heavy flavor production [2]. Again we choose to extract a common factor 
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iCeim^), where Ceirn^) is defined in (2.3). The coefficients of the ^-expansion are denoted 
by i.e. we write 



The ^-expansion of the two-point functions starts at e~^. It turns out that B, 
all i. The first two-point function 



(2.6) 



1 for 



Bi = B{p4- p2,0,m) 



(2.7) 



is real for our kinematics which can be seen by drawing the appropriate Feynman diagram 
for Bi and applying the Landau- Cutkosky rules. The same statement holds true for the 
two-point functions B^ and B4 to be discussed later on. One has 

(2.8) 



ReBt'^-- 


= 1, 


ReSr = 


2-^ln 


Rei?« = 


2ReM°) 


ReSf) = 


2ReSf^ 


Im^i^^ = 


0. 



-t t 



T 



T m2 



1 . , -t 
3 m 



ln3^ + 2Li3(^)+Li3(^) 



T 
m? 



(2.9) 



The second scalar two-point function 

B2 = B{p3+p4,m,m) 



(2.10) 



has both real and imaginary parts. Defining P — {1 — Am?' / sY^"^ and x — {1 — (5) / {1 -\- (5) 
we obtain 



ReS, 



(-1) 



(2.11) 



ReSf 

ReB'i^ 



lmB^2^ 



2 + /31nx, 



s(3 



s(3 



(3 --2 In(^) + - ln2(l^) + 4 ln(l - x) - 2 - x) - 4C(2) - 2Li2(a;) 
(1) 



ReB'i^ ^2ReB\'' + p 



4C(2)ln(^) -^ln3(^)+4C(2)ln(l-a;) -21nxln'(l-x) 

4 



1 ln3(l - x) + 2C(3) - 4Li3(l -x)- 2hh{x) 



ImS^°^-7r/?, ImE^^^^TT/? 



2-ln(^) 



(2.12) 



2Im5^^^ +7r/3 



iln^(^)-2C(2) 
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The remaining three two-point functions have a simple structure: 



B4 = B{p2,m,m) — iC^{m)-; 

e 



e{l-2e) 

iC,{m^) |- + 2 + 4e + + ; (2.13) 

(2.14) 



S5 = S(p3+P4,0,0) = iCe{m^) 



r(2 - 2e) \ w? 



(2.15) 



The results for B^ and B^ in (2.13), (2.14) are not separately listed in the standard format 
which can of course be read off from the relevant expressions (2.13), (2.14). The two- 
point function B^ (2.15) has both real and imaginary parts: 



ReB'f^^ = 1, 
ReSf =2-ln4T, 

Re4'^ = 2ReBf -4C(2) + iln2 4, 

ReSf = 2ReS(^) + 4C(2) In ^ - ^ In'^ - 2C(3); 

6 



(2.16) 



ImSf =7r. 



ImS 



(1) 



2 - In 



(2.17) 



.In2^-2C(2) 



We have done various checks on the above results. First of all, they were double- 
checked, i.e. the results were obtained by two independent calculations. Secondly, they were 
checked numerically by verifying that the original integrals (after Feynman parametrization 
and integrating out the loop momentum, and for Bi and B2 also expanding in e) are equal 
numerically to the final integrals. We have also verified our results by extracting the relevant 
expressions from general formulae given in [9,10]. In particular, our coefficients (2.8) may 
be obtained from Eq. (10) of the first reference of [9] and then using Eq. (2.14) of [10]. 
Our result (2.11), (2.12) can be obtained from Eqs. (2.10) and (2.14) of [10]. Finally, our 
expressions (2.13) — (2.15) can also be obtained from Eqs. (10), (17) and (8), respectively, 
of the first reference of [9] . 

There is one more special case of the two-point integral which is needed for the calcula- 
tion of a self-energy insertion into external massive fermion lines. This integral is used for 
the definition of the fermion mass and wave function renormalization constants in the on- 
shcll renormalization scheme. This specific two-point function is given in the Appendix A 
of this paper. 
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III. THREE-POINT FUNCTIONS 



The one-loop three-point functions are defined by [2] 



/d^Q 1 
(27r)" (g2 - mf) [{q + qif - mi] [[q + qi + ^2)^ - mi] 

The three masses mi, 1712 and 7723 come in various combinations of zero and nonzero masses 
where all nonzero masses are equal to m as before. There are six different types of three- 
point functions Cj {i — 1,2, ...,6) needed for our purposes [2]. They have both real and 
imaginary parts except for C2 and C3 which are real. Again, this can be seen from the 
Feynman diagrams representing C2 and C3 and applying the Landau-Cutkosky rules. Their 
£-expansion is again written in the following universal format: 

Q = t C,{m') {Ici"^) + -^d-'^ + + eCt^ + e'cf^ + 0{e')] , (3.19) 

where the £-expansion now starts at s^"^). Note that the C}"^^ are purely real. 

It turns out that the 0{e'^) results for the three-point functions can no longer be pre- 
sented in terms of classical polylogarithms but require a new class of functions given by 
the one-fold integral representations defined below. To write down our results in a short 
and convenient form, we introduce the following functions: 

T I \ ln(Qii + aiy) ln{a2 + a2y) Inja^ + a^y) 

JQ 014 + y 

and 

T ( \ /-^ , ln(ai-Haii/)Li2(a2 + tt3|/) .oo1^ 
L(^i(q;i,q;2,q;3,Q;4) = / dy ■ . (3.21) 

Jo Oi^ + y 

Here the ai {i = 1,2,3) take values ±1 and the a^'s are either integers {1,0,-1} or else 
kinematical variables. The above L-f unctions arise naturally in our calculational framework 
^. They can all be expressed in terms of so-called multiple polylogarithms of maximum 
weight four [14] (see Appendix B for details). However, we choose to write our results 
in terms of the above single- and triple-index L-functions for several reasons. The results 
look simpler, e.g. they can be expressed as one-fold integrals of products of logarithms and 
dilogarithms, and are shorter. We have also found that the L-f unctions are much easier to 
evaluate numerically than the corresponding multiple polylogarithms (see [15] for relevant 
details). 

There exist simple algebraic relations between these L-functions based on either sym- 
metry relations regarding permutations of indices and change of integration variables or 
on relations based on integration-by-parts techniques. We describe them in Appendix C. 



^As A. Davydychev informs us the functions analogous to our triple- index functions L^^„^„^ also arise 
in the approach of [10] when one analytically continues their Eq. (3.2) for the order terms. 
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In particular this means that our results on the three- and four-point functions can all be 

written in terms of the L_^^ and variants of the triple-index Lo-^ 0-20-3 functions in 

Eq. (3.20), and of the L+ variant of the single-index Lg-i function of Eq. (3.21). 
We start with the three-point function Ci defined by 



Ci = C(p4,P3,0,m,0). 



One obtains 



(-1) 



(3.22) 



ReCl 



(0) 



ReCi 



(1) 



ReCl 



(2) 



1 



2s(3 



[ln'x + 4Li2(-x)+2C(2)], 



1 

'7^ 



- In^ — + 2 In — ln(l -x)lnx + ln(l - x) In^ x - 4C(2) In — 



6 m 



+ 21n— Li2(x) +5C(3) -4Li3(l -x) + 2Li3(a;) + 2Li3(l - a;' 



- SLia 
1 



■1 + x' 



7,4s 11,3 

— In^ — In"* 

48 



1 s / 

Inx In — - (In (1 — a;) -|- 61n(l — a;) In 

4 777,2 V ^ 



X 



+ 5 In^ x) + In 
1 



7 



11 



^•^(1 -x) In (1 - x) \nx H x 

V ^ 2 ^ ^ 24 



^1 5 55 

— - ln^(l — x)lnx ln^(l — x) In^ x -\ — ln(l — x) In^ x -\ In^ x 

3 4 2 48 



+ 



31n^ 



10 In 



lnx + 161n(l -x)lnx + 31n^x C(2) 



In 



11 Inx) C(3) + I C(4) - 2Ul(-x) + Lhix) In^ 



+ 21n(l -a;)lna; + 71n2a;-21n— (ln(l -x) - Inx) - 10C(2) 



13 , 



+ Li2(-a;) In^ — - 2 In — (ln(l - x) - Inx) + 21n(l - x) Inx 

V 4 

+ — ln^x-6C(2)^ +4Li3(l-x) fin ^ + 2 ln(l - x) - Inx 
4 / \ m"* 

/21 9 29 \ / s 

- Li3(-x) — In — + 201n(l - x) ^ Inx - 4Li3(x) In — 21n(l - x 

\ 2 2 / \ 

+ \ ln(x)) - Li3 (^) (3 In + 4 ln(l - x) - 2 In x) - 2Li3 (l - x^) x 



In — 21n(l - x) - Inx ) - 4Li4(l - x) - 4Li4 



-X 



1 -X 



+ 12Li4(-x) 



+ 28Li4(x) + 17Li4 (^) - 17Li4 (3^) + 4Li4 (l - x^) + 4Li4 
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- L_++ (l, O) + L_++ (l, + 2 L_++ (l, X, x'"^ O) 

- 2 (^1, x, a;""^, a;~^^ — 2 (l, a;""^, a;^ + L_++(l, x, 0) 

- L_++ (l, X, X, x"^) + 2 L_++ (l + x-\ 0, 0, -l) - 2 + x, 0, 0, -1) 

- 2 L_++ (l + x'S 0, 0, -1 - x) + 2 L_++ (l + 0, 0, -1 - x"^) 
+ 2 L+ (O, 0, x) - 2 L+ (O, 0, -X, x'^) - 4 L+(0, 0, -x, x) 

- 4 L+ (O, -l) + 4 L+ (O, -1 - x) + 4 L+(0, -x, -1) 
+ 4L+ (O, -a;"\x"\ -1 - x"^) - 4L+(0, -a;,a;, -1 - x) 

- 4 L+ (O, -X, X, -1 - x-i) - 4 L+ (o, -l) 
+ 4 L+ 0, , -1 - a; + 4 L+ 0, , -1 - x'^ 



V 1 — a;'l — x^' 7 V'l — x'l — x^ 

+ 2 L+ (x-\ 0, -x-\ x) + 2 L+ (x-\ 0, -x, x) - 3 L+++ (o, 0, x'S x) 

- (^0, 0, X, x~"^) + - (^0, x~-^, x'^, x) + 3 L+++(0, x~\ x, x) 



2 

~t~ — I — [_ ^0, X, X , X ) 2 — ' — ^ — ^ — ' — ^ ' ' 



TT 



ImC}~^^ = 0, ImCf^ = ^Inx, (3.23) 



2s/3 

TT 



2 In — In X - 4 ln(l - x) In x + In^ x + 4C(2) - 4Li2(x) 



6sp 



3 In — - In X — 12 In — r- ln(l — x) Inx + 31n — -In x — 61n(l — x)ln x 



+ In^ X - 12 In 47Li2(x) - 24Li3(l - x) - 12Li3(x) + 12C(2) In ^ + 12C(3) 

We have not been able to derive the corresponding result from the known general hyperge- 
ometric function that represents the above integral in [9]. On the other hand, for a general 
three-point function, an expression for the order £-terms was obtained in [16] in terms of 
simple polylogarithms up to Lis. However, we believe that the result Eq. (5.21) in [16] is 
not applicable to our case as one faces singularities resulting from vanishing denominators 
in the arguments of the relevant logarithms and polylogarithms. We have checked our final 
result numerically against the original two-fold and one-fold Feynman parameter integrals 
(after ^-expanding the corresponding integrand). This was done term by term for coef- 
ficients at the corresponding orders in e. Although the result for the £^-coefRcient looks 
lengthy, our final analytic results (3.22), (3.23) for the three-point function Ci integrate out 
numerically very fast (in fraction of a second on a desktop computer for a chosen numerical 
point) and without any problems. In comparison, the numerical integration of the one-fold 
integral by Mathematica took eight times longer, and that of the two-fold integral even 200 
times longer to evaluate. In addition, because of various branch cuts, the one- and two-fold 
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integrals would only allow integrations in the complex plane of kinematical variables, while 
for the physical region they have severe problems. 
The integral C2 is real and finite: 

C2 = C(p4, -p2,0,m,m); 



ReCf) 
ReC« 



ReC?^ 



= ReC^"^^ = 0, 



(3.24) 



1 

I 
1 

M 
1 

m 



C(2) - Li2(^ 



-t 



-t. 



In' — + 61n — Li2(— ) + 9C(3) - 6Li3(-) - 9Li3(— ) 

t 

-t. 



-5 In^ 4 + 12C(2) In^ 4 + 12C(2) In^ 4 



121n2^Li2(4) 



-24 In 4 In 4-Li2(4) " 24C(3) In 4 + 24C(3) In ^ 

+241n4Li3(^) + 241n4Li3(T) + 241n4Li3(4) 
-T T 

+481n— Li3( — ) + 192C(4) - 24Li4(— ) 



ImC^^') 



+12L_++(1, 
0. 



2 2 



-t 
t 



^,0)-12L_++(-, 0,0,-1) 



(3.25) 



This result was checked numerically against the original double parametric representation 
(obtained after doing Feynman parametrization) of this integral expanded in powers of 
e. We could not obtain similar expressions from known general results for this integral, 
as the £-expansion of the relevant hypergeometric function is problematic. In addition, it 
turns out that the general result for the order £-terms for the massive three-point function 
of [16] does not allow for a straightforward extraction of the corresponding expression for 
this particular case. More exactly, the equation Q:j,{y) = originating from the table in 
[16] (on page 608), does not have solutions for the relevant kinematics. In this sense, our 
expressions for the coefficients of the e- and £^-terms for C2 represent a new result. 
The integration of the function C3 defined by 

C3 = C{-p2,P4,0,0,m) 

requires the construction of a subtraction term since an e-expansion of the relevant inte- 
grand does not straightforwardly lead to the desired ^-expansion of the integral. This is 
best illustrated in a simple example which nevertheless captures the essential idea of the 
subtraction method. Consider the integral 



/ dxx ^^^f(x,e) 
Jo 



(3.26) 
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where /(x, s) is an integrable function in the interval [0, 1] and has derivatives in e. For the 
sake of the argument take f{x, e) to have a Laurent series expansion starting at the zeroth 
order in e, i.e. f{x,e) = f^^\x) + ef^^\x) + It is clear that expanding the integrand 



X 



X \ X X 



x)]+... 



does not in general render the integral integrable. However, if we write 

£ dxx-^+'f{x, s) = dx {x-^^')^ f{x, s) + /(O, e) dxx- 



l+e 



(3.27) 



the terms on the r.h.s. of (3.27) are now integrable. In (3.27) we have introduced a "plus" 
prescription 

£ dx ^ fix, e) = £ dxx-'+%f(x, s) - /(O, £)) 

not unlike the "plus" prescription usually introduced when discussing parton splitting func- 
tions. The £-expansion of the integral (3.26) can now be obtained since the first integrand 
on the r.h.s. of (3.27) can be expanded in e and then be integrated term by term whereas 
the second integral can be computed in closed form. The task is then to find the appropri- 
ate subtraction terms for the integrals encountered in our calculation. This is required for 
the three-point function C3 and the three four-point functions to be discussed in the next 
section. 

As exemplified above we derive the subtraction terms by substituting the value of the 
integration variable (usually the lower or upper limit of integration) at which the given 
integrand diverges into the nonsingular part of the singular integrand. Adding and sub- 
tracting the subtraction term does all the job: e.g. the subtraction term contains all the 
poles in a given Feynman parameter but can be easily integrated due to its simpler analytic 
structure, while the rest of the integrand is now finite with respect to the same parameter 
and can therefore be integrated as well. When dealing with such a finite but complicated 
integration we often make use of the integration-by-parts method to evaluate and simplify 
our expressions. 

Applying the subtraction method to the evaluation of the three-point function C3 one 
obtains: 



Re a 



(-2) 



1 

2t' 



ReCt^ 



ReC« = 



t m'^ 

^ln3 4 + 2Li3(7) + Li3(^; 
3 m'^ t 



ReCi'^ 



-t 



T 



(3.28) 



3 — 



In^ 



-t 



-In' 



-t 



In 



-T 



1 

M 

T w? T 

-3Li4(^) + 6Li4(— ) - 6Li4(- 



2-3C(2)ln- , 



V6C(3)ln4-6C(4) 



-t 



t 



0. 



(3.29) 
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Note that one can obtain corresponding expressions in terms of generalized Nielsen polylog- 
arithms from Eq. (27) of [9]. The corresponding hypergeometric function of three variables 
$1 can be reduced to a hypergeometric function 2F1 of one variable and one can then use 
Eq. (2.14) of [10] to get the relevant ^-expansion. We have verified agreement with [9] 
analytically up to 0{6). The agreement for the £^-terms was verified numerically. 
The three-point function C4 has a closed form solution: 



s + 



s r(i-2£) 



(3.30) 



which is straightforward to obtain. For the £-expansion of the real and imaginary parts of 
(3.30) we get: 



Re a 



(-2) 



(1) 



Re a 



ReCP 



Im Ci"^^ 



ImCf 
ImCP 



ReC, 

4C(2)ln 
1 

24 



(-1) 



^1 ' 



Reef 4 



In' ^ - 8C(2) 



1 



In^ 



s 



- 2C(3) 



In' - 2C(2) In^ ^ + 2C(3) In — + 9C(4) 



s 



s 
n 

71 



(3.31) 



(3.32) 



In^ 



4C(2) 
1 



2C(2)ln4--ln3 4-2C(3) 



For the fifth three-point integral C5 defined by 

C5 = C{-p2, -Pi,m, m, m) 



we first obtain a one-fold integral representation similar to Eq. (3.13) of [10]. As before, 
the main difficulty is the derivation of the coefficient for the £^-term. The corresponding 
coefficient has a complicated singularity structure as well as two branch points on its 
integration path. Therefore, in order to analytically separate the real and imaginary parts 
for our final result, we have divided the integration regions for the relevant terms into 
three parts. After analytical integration these terms are free of numerical instabilities and 
converge very fast. One obtains 



ReC^^^ 


= ReCi; 


ReCf 


1 


In 




2^ 




ReC^^^ 


1 


T 




Ys 


.3 



0, 



(3.33) 



6C(2)] , 

- 8C(2) In a; + 4 In xLis (x) - 6C(3) - 8Li3(a;) 



12 



1 



9 11 

41na;ln^(l - x) ln^xln^(l - x) ln^xln(l - x) ^ In^ x 

2 3 12 



-C(2) (31n2x-101nxln(l-a;) + 61n2(l-a;)) xU2{-x) 

+2 (in^x - 61nxln(l - x) + 31n^(l - x) - 6C(2)) Li2(x) + 10C(3) \nx 
-6C(3)ln(l -x) +2(lnx + 31n(l -x))Li3(x) + 21nxLi3(-x) 

-12 (Inx - ln(l - x)) Li3(l - a:) + y C(4) - (l, 0, 0, 

1 — X 



+L_++( 1,0,0, 
-L 



0,0, 



'+++ 



1 — xl — a; \ 
0,— ,— ,-1) + ^+++ 



X 



X 



X 



0, 



+++ 



1 — X 1 — X 1 



0,0, 



1-x 1 



X X 



X 



X X 



ima 



(-1) 



Ima 



(1) 



^ 6s 



0, ImC^°^ = -Inx, 
s 

TT 



(3.34) 



ln^a;-4C(2)+4Li2(a;) 



-121nxln2(l - x) + In^x - 12(lnx - 21n(l - x))(C(2) - Li2(x)) 



+12C(3) - 12Li3(x) - 24Li3(l - x)] . 

Explicit result for this integral was given very recently in Eq. (4.4) of [11]. We have checked 
agreement with [11] analytically up to 0{e). The agreement for the £^-terms was verified 
numerically. 

Finally, we write down real and imaginary parts for the last required three-point function 
Cq defined by 

Cg = C{p3,pi,m, 0,m). 



One has 



Re a 



(-2) 



0, Red = ^Inx, 

sp 



(3.35) 



ReC, 



(0) 



1 



(1) 



ReC, 



ReCP 



2sP 
1 



-41na;ln(l -x)+hi^x- 8C(2) - 4Li2(x) 



-6 In^ X ln(l -x)+ln^x- 24C(2) Ina; + 72C(2) ln(l - x) + 12({3) 
-12Li3(x) -24Li3(l-x)], 



1 



24s/3 



-161nxln^(l -x) + 241n^xln^(l -x) - 81n^ j;ln(l - x) + In^ x 



+41n^(l -x) + 192C(2) lnxln(l - x) - 48C(2) In^ x - 240C(2) ln^{l - x) 
-48C(3) Inx + 120C(4) + 48Li4(x) + 96Li4(-— ^) + 96Li4(l - x) ; 



TT . ^(0) TT 



— , ImCr = ^[lnx-21n(l-x)], 



(3.36) 
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ImC^^^ = — [-121nxln(l-x) + 31n^x + 121n2(l-x)-12C(2)], 



g _ — ^121nxln^(l -x) - 61n^xln(l -x) +ln'^x -81n'^(l -x) - 12C(2)lnx 
+24C(2)ln(l-x)-12C(3)]. 

Corresponding results for Ce may be obtained from Eqs. (3.5), (3.7), (2.10) and (2.14) of 
[10]. We have done an order by order numerical comparisons for the coefficients of the 
£- and £^-terms, while other terms can be easily compared analytically. We obtain exact 
agreement. 

We mention that we have checked all our analytical results for the three-point functions 
against numerical results provided to us by M.M. Weber [17] (see also [18]). We have found 
agreement. 



IV. FOUR-POINT FUNCTIONS 

The scalar four-point one-loop integrals with one, two or three heavy quarks running in 
the loop are the most difficult to evaluate. The one-loop four-point functions are defined 

by [2] 

Diqi, q2, qs, mi, m2, m^, 7714) = (4.37) 

2e f 



J (2^ 



(27r)" (^2 _ ^2) [(^ ^ _ ^2] [(^ ^2] [(^ + + + _ ^2] ■ 

As before, the -\-i5 terms in the denominators have not been written out. Again, there is 
only one internal mass scale for our purposes. 

For heavy flavor production one needs three different types of four-point functions 

Di (i = 1, 2, 3) which arc expanded as 

A = i C.(m^) [^d\-^^ + -^d\-'^ + Df + + e^of^ + 0{e^)^ . (4.38) 

Again the coefficient of the most singular part of the four-point functions is purely real, i.e. 
Im^H^ =0. 

Before we give our results for the four-point functions it is necessary to discuss some 
general technical fcatTircs. After applying Feynman parametrization, we are left with a 
three-fold parametric integral for the four-point functions: 

D = iCe{m^)il + eMm^y^ C dxidx2dx^ xlx2K-^~', (4.39) 

Jo 

with the kernel K given by 

m^K = —abql — ac{qi + ^2)^ - ad{qi + q2 + q^Y — bcql 

—bd{q2 + qzf — cdql + am\ -\- hm\ -\- cm\ -\- drril — i5, (4.40) 
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where the (i = 1, . . . , 4) can be either m or 0. The set of parameters {a, b, c, d} above 
corresponds to an arbitrary choice from the set of original parameters {X1X2X3, a;ia;2(l ~ 
X3), Xi(l — X2), 1 — Xi}. For each particular four-point function we make a judicious choice 
of Feynman parameters which enables us to get the most convenient kernel for the subse- 
quent integrations. 




^2 ^4 
FIG. 2. Massive box Di with three massive propagators. 

First we consider the four-point function Di with three massive propagators shown in 
Fig. 2 which is defined by 

Di = D{p4, -p2, -pi, 0, m, m, m). 

Substitution of the corresponding values of momenta and masses for into the expression 
for our kernel (4.40) gives 



K — act — hds + (1 — a)^ — 
where we have introduced positive valued dimensionless variables 



t 



(4.41) 



(4.42) 



The kinematical conditions s>4,f>l,s>f constrain the allowable region of phase 
space for the present physical 2 —>■ 2 process. Our choice for the parameters {a, b, c, d} is 
{1 — xi, X1X2XS, 2:1X2(1 — X3), xi{l — X2)}- For Di, the integration of the corresponding 
integrand over 0:3 results in two terms: 



X1X2 



JJD, 

X1X2 



X, 



Xi + t{l — Xi)x2) 



-1-6 



{l + e)[sxi{l-X2)+i{l-xi)y 
{l + e)[sxi{l-X2)+i{l-xi)] ' 



(4.43) 
(4.44) 



which then have to be integrated over the remaining parameters xi and X2- Eqs. (4.43) and 
(4.44) correspond to the indefinite integral (or primitive) evaluated at the upper and lower 
limit of X3, respectively. The term I^^^.^ in (4.43) does not change sign on the integration 
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path, e.g. does not have a branch cut in the interval [0, 1] for both variables xi and X2- 
Consequently, it does not give an imaginary contribution and it is thus safe to drop the 
iS shift in IxiX2- Furthermore, since I^^xi does not have poles in e, we expand it up to 
£^ and straightforwardly integrate over the second variable X2 to obtain Concerning 
the second term in (4.44), one can see that there is a branch cut for the variable X2 

in its numerator as well as a divergence due to the factor Xi^~'^^ at the lower limit of the 
integration Xi = (wc have dropped the iS shift in the denominator as it does not affect 
our further calculation). At this point we introduce a subtraction term for 11x^^2 
simplest possible way: we set xi = in //^^j everywhere except for the divergent term 
^-i-2£ rpj^-g pgg^i^g ^jjg following subtraction term: 



_ Xr^-^" [1 - g.T2(l - .X 2) - i5]-^-' 

'I + e)t 



IIxX = M ' (4-45) 



which, in the framework of the dimensional regularization scheme, integrates over xi to 

2e{l + e)i ■ ^ ^ ^ 

Then we expand the above expression up to e"^ and reexpress the argument of subsequent 
logarithms as 

{x2- xP){x2-l+ xf^) 
xf(l-4°^) 

with 

1 + Jl - 4/s 1 + i3 

4'^ = 4 +tS = ^ + z5 (4.48) 

A final integration of the subsequent series can be done analytically in the complex plane 
and its result is expressed in terms of logarithms and classical polylogarithms up to Li^. 
Analytic continuation of the result for 5 — > is then straightforward. 
Lastly, we calculate the finite difference 

{l + e)i[sXi{l-X2) + i{l-Xi)] 

by again expanding the difference up to £^ and using (4.47) for the arguments of the 
logarithms. Then we first integrate over the variable X2, leading to a reduction of the 
integrand to simple fractions w.r.t. X2. In this way we avoid spurious poles in the remaining 
integral which would otherwise arise in case of integration over xi first. 

To complete the derivation of the first four-point (box) integral, we combine the two 
terms 

Ix,' + A//i^^ 

and perform the last integration over the variable xi. 
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At this point we would like to comment on some technical details of our calculation 
which are used throughout this work. For instance, the integrand for the last integration 
contains expressions such as 



f{xi) ■ Li2,3 



aixf + a2Xi + as 
a4xl + a^Xi + Qq^ 



(4.49) 



where f{xi) is a rational function or a product of a rational function and a logarithm. 
Using recursively the method of integration-by-parts as much as necessary we render their 
arguments to be linear functions of xi. In addition, in the case of Lis, we can reduce the 
weight of Lis by one. At the same time, the sources of imaginary contributions are trans- 
ferred into logarithms (or remain in Li2's and Li3's with arguments that are independent 
of the integration variable) . Finally, performing the last integration and adding up all the 
relevant contributions we arrive at the result for the box integral with three massive lines, 
containing polylogs up to Li4 and the single- and triple-index L-functions introduced in 
Eqs. (3.20) and (3.21). 

In order to keep our results at reasonable length we introduce the abbreviations 





Z3 
Z5 




Is 




h 


and obtain: 




ReD\'^ = 


0, 

- 21nx 



Z3 = {s + 2t + sP)/2, 
t + tp)/2, 



= {2im? 
ln/3. 



k 

lz3 ~- 



In 



-t 



Z4 = (s + 2t - 

Zq = {2rn? - 
-T 



li 



In- 



In 



2 ■ 



= In 



-za 



- s(3)/2, 
t - 1(3)/2, 

Ij. = In X, 



(4.50) 



m 



ReD[~^^ ^\nx/{stp), 
21nxln(-t/3/m^) +2Li2(x) -2Li2(-x) +3C(2)] /{st(3), 



(4.51) 
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k {ll + II3 + 21t (Ix + lz3 - Iza) + + Ix (2 lz3 - 4 k)) + Ix I, 

Is (it Ix + 11/4: + Ix {lz3 + Iza) - if) + (-/. + m -lx+ 6h) C(2) 

5C(3) + 4Li2(x) Zt-2Li2 



rn? X 



2U2 



2w? 



k - 2Us{-x) - 2Us 



-X 



-2Li3 



+ 



-3.f)..U3(i^)..U3(^ 



2^5 



8Li, 



13' 



2P 



2Li, 



' 2z, ' 
m2(l + /3^ 



+ 2Li3 



2zq 
-i(l + /3). 
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1 



ifi _ 1^ 73 7 _ Z /2 72 2 3 5^ _ 11 

"64 24 4^.^* + 3^3^t + ^2 24 



IIIt- - ll It h- 



ib 4 2 

4 Is If Ix ~ h ~ o ^tIx ''r — h h It^x ~ h + q h ^t^x J' h 



13 
T 



3 
2 



3 



3 3 17 3 7lt 29 
o ^-+24^^^- ■ 64 ■ 48 

— It It lz3 + 2 + 2 ^^3 ~ ^a; ^28 + h k Ix hs + 2 hs + 



11 



1 

8 



1 



^ lz3 ~ g ^z3 ~ 4 ^23 ~'~ 2 ^* ^^"^ ~'~ ^^^^ — 5 It It 1% 7 ^o; l^-i + 

27 .0 , . 1 



^ I's I'x I'z3 



2 2I22 3 61 3 2 - -1 

k Ix lz3 — '^ItIx lz3 ~ Q^x ^23 + h hs 4g ^24 + — h k hi + ^^^^t hi 

— If lz4 + -I^It hi — -^hk It hi + 2 hi + h ^T hi ~ ^h^T hi + 



9 11 ^ ^ ^ ^ 

- ll h hi + -lslt h hi + - ll h hi --hh h hi --klT h hi +-lTh hi - 



16 
11 



3 

2 



3 
2 



5 
2 



1 

2' 



9 



-j^g h Ix hi ~ ^Itl^ hi — '^^tIx hi ~ Yg ~ ^ ~ ^ ^'^ ^* ^'^'^ ~ 



I'S H I'z3 I'zi 
1 



It h3 hi ~ 3/s It hs hi ~ It hs hi — h hs hi — 3/r Ix lz3 hi — lz3 hi + 

5 



77 



13 



1 



2 It It lid + T ^24 TT h Ix lid ~ '^hh lid + — It h ltd T7^ It lid ^ h lz3 ltd + 



33 



5 



7 



h lz3 hi + o ^* ^23 hi + h lz3 hi -.r^s hi o h h hi "I" 3 Z 4 h It hi + 

2 lb 2 4 
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2 ,2 



8 



4 
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h lz3 lid + 3 h3 lid + 2 ^23 ~ ^23 ^24 ' ' '3 I ' '3 ; ;3 



— Is lid + 2 ^24 ~ g It ltd + 



12 



49 



13 



9 



X •'24 



2^23 hi ~ —l-yA — ITT h Ifi — — h Itlfi — 2 If, h Ift — — h liH r h It I 



24 



^24 



24 



1 1/3 



13 1 

3 It 1/3 — 4: h ItIp + — h hh — '^ h k hh — ^ k hip — h It h h ~ 

o 2 

2 9 \ 7 17 9 

3 It It h h + lTlxh - h lxh+ 2 ^* 4 24 4 ^'^ 



4 is ^23 h + 4:lt hs h + '^IsIt lz3 h + (ilt It lz3 h ~ 2 ^-^3 h + 4ltlx lz3 h + 

5 11 

4^r ^23 h + ^23 — h lz3 h ~ 3it h^ h ~ h hs h + ^24 h + ''^s ^^4 h + 
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1/3 + ^Is W IzA 1(3 + '2 k It IzA + 2 + 3 Ix IzA 1(3 ~ hi ^(3 ~ 

- ll IzA — ha hi ~ ha hi Ip — 4:It hz h4 1 13 + h hs h^ + 2 hi 1(3 ~ 
5 Is h ~ h ~ '^h 1(3 + ha h + 2^z4 h ~ g^s ^0 ~ l^h hip — ~ 

-hlrl'p — ^ h ItI'p + ^ h hl'p — 2 ^* hl'p — 2 hl'p — -^Ix^'p '''' 2 ^"p 

2 2-'^22252 2 2 

3 h ha lf3 + ha l(3~ i^h ha 1(3 ~ ha h ^ + 3 /t h4 Ifj + 21t hi 1(3 + 

— Ix hi I'p ~ 2 hi ^"p ~ 2 ^'p ~ '^hlp ~ '^hl^ ~ It ~ g h Ip + ha + 
Ib /19,o „,o 11. 



hJp "g' ( "g"^^ ~ "I" ~ + ^hh + hha + -^hha — ^la ~ ^hh4 + ^hhi— 



3 ha hi -7ll^ + 2lt {31t + 7 ha + 5 - 12 h) + Qhh - l^hk - 2 ha h + 
14 U //3 - 13 /| - J (64 k-8lT + 35h + 18 ha - 44 U + 32 /^)) ((2) + 

(-2 - 4 - Zt + + 7 + ^.4) C(3) - yC(4) - 2 Li^ f — j + 

2ll + l^ + lTh + f + 2lTha + Shha + k (-4Zt -Sh + 21,^) + 

ik -It- - ha - h^ + 5 L. h^ + Ha - A L - 6 C(2)] + U2 (— 1 x 



2 ^.4) h4 + 4 - 4 /c. //3 - 6 C(2)) + Li2 (-^^ 

f -2Li2 ( — ) - f + + + 4^' - f (6^* + + 3/. - 4/,3 - 4/^4) + 

2 h h4 - IIa -lt{h-2lT + 2 ha + 2 h4) - 2 - 6 C(2)) + 

Li2 {^^^^^) (-2Li2 (J) + \ [l] - I2ll - 2ll - 6hh + l2lTha+ 

8 h ha -4.lxh4-4^ll4 + h {-4^k + QlT + 2h-S ha + 4.1,^-41^) + 
4 (h + 4 + 2 - 2 Ip) + 121x1(3 + 8 h4 h-4 if) + 12 C(2)) + 

1 Li2 (^^^ (-II + IQlth - h (4/t + h- 4/.4) - 2ls {21t + h- 2/.4)) + \ 



X 



Li2( — ) (3/,' + All + 2/t + 2lTh + 7ll + Slrha - 2/^(8/* + Sh -h + Aha - 6/m) 

+4/t^z4 + Alxh4 — 12/^4 + ^kiha + h4 — 1(3) + ^1x1(3 — ^hal(3 + 8^x4^/3 — 4i|) + 
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^ Li2 (^^) (9 - 16 + (-4 It + 1. + 12 /,3 + 8 /,4 - 8 /^) - 

2/,(4/t - 2/t + 5/^ + 6/,3 + 4/,4 - 4/;?) - 8/t(2/T + k - 21,3 - 21,4 + 21 p)) + 

\U2{^) (211 - 2lTk - 21x1,3 + ^Uz3 + - 2lt{k + 21,3) + 21x1,4 + ?>IJ,4+ 

,2 - - - ^- - 1- 



^lzzlz4 - 1,4 - 2Uf3 - Al,3lp + ls{2lt - 3h - 71,3 " lz4 + 21,3)) + -^U2{-x) X 

(5/2 - 16/2 + ll- AItI,3 + QUz3 - GII3 - ^ItU + 21J,4 - Sl,3l,4 + 2/^4 - h+ 
8l,3h + 2ls(-Glt + 21t - 21^ + hs - Sl,4 + 2lp) - Alt{l^ - 61,3 - 4Z.4 + 4Z^)) + 

Li3 (I) {-It + 1,4) + 2(1.13 (5)+ Lis (g)) {2k + 21t - 1,3 - U) + 

2U3 (^^) {2h-2lt-2lT- 1,3 - U) + 

Lis (^-Y/.-8Zt + 4ZT- | + 8Z,3 + 9Z,4-8Z^j + 

^ Li3 (^^^ ( -5 - 4 - 2 + + 6 + 4 ^,4 - 4 i/j) + 

2 Li3 (j-^ {-Is + lt + lT + lx- h) + Lis (y) - 2 /t - 3 + 2 1,3) + 

ilia (^^) (Z. + 8Zt - 81t -k- ^ha - ^U) + U3 (^y ) {h + + 21t - - 

21,3 - 4U) + ^Li3 (I) {Is + 21t + k- 21,4) + 2 (^Li, (j^) - 

-3 ( V) -K^) -K^) -3 Q) - 

/t - lz3 - IzA + - 2Li3 (J) - Lis (-) lz3 + (Lis (-) + 2 Li3(-a;)) 
(2Z, + 2lt + It- 21,3 - 21,4 + 21 p) + 2Lis(x) (4/, + 3/* + /r + - 2/,s - 6/,4 + 
3Z^) + 2Li4(x) - Li4 (y ) - 4Li4 (I) + 4Li4 (|) - Li, (|) + Li, (|) + 

(^) (4^) . .LU (Zl±i!) . (^) . 

1 ( s{l-(5)\ 1 ( rri^ s{l+(5)\ 

-L_++ II, — , — , — ^ I - -L_++ II, — , — , I + 



X 
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^3 ' ^4 ^ 
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/ s{l-(3) s{l-(3) \ _ 1 / s{l-(5) .(1-/3) .(l + /3) \ 
V ' -2-^4 ' -2-24 ' / 2 1^ ' -2^4 ' -2^4 ' -2^3 y 

0. 0. -1) - iL_,, (1. 0. 0. -1) + 2L^,. . »■ -0 - 

It' °' "i^' ~VJ " It' °' "i^' " It' "i^' "i^' 

,i ,3 /t i\ ft t \ 

3L_++(-, 0, 0, -1) - -, 0, 0, 2L_++ -, 0, , -1 + 

Z4 z \Z4 zs/ \Z4 Z3 y 

ft t t\ ^ f t ^ t t\ 
L_++ -,0, , +3L_++ -,0, , + 

\Z4 Z3 Z3J \Z4 Zs Z4J 

1 fttt t\ f sjl + P) m'\ 

2^-++ U' "i^' "i^' -f) + 2^-++ V^^' ^' ' " 

2^3 -r' -2^4 ; 2^3 -r' -2^3 

2^^ ^' ^' ) + Iv^^' ' ^ J 

1 f s{l + f3) s{l-f3) s{l-f3) m'\ ( z, z,z, \ 

2^-++ ^2ir' ^2ir' ^ j " 1 ^' ~7 

4r /"n ^3(1-/3) ^3^4 /m^ ^3(1-/3) ^3^4 s{l-ff) \ 

V -2t/3 ' 1^' + l,^' -2t/3 ' ^2ir j 

2L, (-1, 0, ^, -1) + 2L, (-1,0, f , -1) + 2L, f-^, ^, 

V 2^3 ^ ^4/ \ Z^ t Z4/ \ Z3 S/J StjJ I ^ 

2L+++ ^0, — , — — ^ J + 2L+++ ^0, — , — — ^ J + 
/m^ g(l + /3) ^ 1 /m^ g(l - /3) s{l + P) \ 

V 2^3 / ^ V 2^3 2^3 Z4J 



+ 
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'+++ 



' s{l-(5) s{l-(5) 
—2z4 —2zi 



,0, 



m 



-T . 



lmD\ 



(-1) 



^ stp 



7r/{stp), 

TT 



ImD[°^ = -27rln(-i/3/m^)/(si/5), 



(4.52) 



3 2 2^2 2 

^ + 2 ~ ^ ^a; ~ ^xhs + ^a;^z4 ~ ^24 + h {k + ^2:/2 + lz3 + 



{Ix + 2k3-Alp) + 2ll- 2C(2) - 2Li2 



m 



2Uo 



2m? 



TT 



-(7Zf-Z^-3Z,2(3Z, + 4Z,3 + 4U-6Z^)+6Z^(-2Z,3 + 3Z;,) + 



12 (/,\ - 2 (-/23 + 44) //3 + 2 (-43 + Ifi)) - 
8lp {6l^t-3ll^ + 2ll + 6k{-ks + h))-Sls (8/2 -5/2 -4/^4- 
8 /* {Ix + lz3 + 8 lz3 h + ^ h + ^lx (-2 /.a + 21,4 + h))) + 
4 ik + lx + h) C(2) + 2 Li2(x) k + 2 Li2 k + 

2 Li2 ( ^^~/J"2^^ ) ih -lx + 2l^) + 2 Li2 j (/. + k + 2l^)- 

2L.3(.)-4L.3(f)+2L,(|)-4L.3(f)-2L,g 




FIG. 3. Massive box D2 with one massive propagator. 

Next we turn to the second four-point function D2 with one massive propagator shown in 
Fig. 3 which is defined by 

D2 = D{-p2,P4,P3,0,0,m, 0). 

Wc substitute the appropriate values of momenta and masses for the D2 integral into the 
general kernel expression (4.40) and obtain 



K — act — bds + (? — iS. 



(4.53) 
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In order to simplify the first integration over the Feynman parameter we choose {a, b, c, d} 
as {xiX2{l — Xs), X1X2X3, Xi{l — X2), 1 — a;i}. After a;3-integration we write the result for 
the integrand as a sum of two parts /^^^ + ^^xiX2 '■ 



Xi 



xj{l-X2)il + {i-l)x2)-i5] 



{1 + e)[S{l - xi) + txi{l - X2)] 
jjD2 ^ a;i [-SX1X2 + xjil + (g ^ 2)3:2 + xl) - i(5]~^~^ 

(1 +£)[§(! -Xi)+tei(l-X2)] ' ^ ■ ^ 

where again the two terms derive from the indefinite integral (or primitive) evaluated at 
the upper and lower boundary of X3. First consider the integration of /^^^j- One notes 
that its numerator is not negative on the integration path since t > 1, which implies there 
is no imaginary contribution coming from Z^^^. Therefore, we omit the i5 term for the 
remaining integration. After integration over xi we arrive at 



^j,, ^ [(1 - X2) (1 + {i - 1) X2)] '''' 2Fi(l, -2£, l-2e,l-A) 



X2 



2s (! + £)£ 



where we have defined A = t {1 — X2) / s and 2-^1 is a hypergeometric function. The above 
expression is singular at the upper integration limit X2 — 1. In order to regularize this 
singularity we have to find a suitable subtraction term. 

First note that the ^-expansion of the hypergeometric function reads 

Fi(l,-2£, 1 - 2e,l- A) ^l + 2e \n A - Ae'^ U2il - A) - 8£^Li3(l - A) 

-16e^Lu{l- A) + 0{e^). (4.57) 

To obtain a suitable subtraction term we substitute (4.57) into (4.56) and replace X2 by 1 
everywhere in except for terms that diverge. Therefore, our subtraction term can be 
defined as 



.D2,s ^ (1 - X2)-'-H~'-^ [1 + 2£ In A - 4eX{2) - Se^ ({3) - IGeXm 



2s{l + £)e 



(4.58) 



Our subtraction term is simple enough to be integrated analytically over X2 giving the 
result 

^ _ (-3 - 2£ In I + 4e^ ({2) + 8e' ((3) + 16s' C(4)) 

2{l + s)s^s ^ ' 

which can readily be expanded in e. 

Next we turn to the remaining finite integral I^^ — I^^'^ ■ We expand I^^ — I^^'^ up to 
second order in e and integrate over X2 after the expansion. 
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Next consider the second integrand //J^^^ (4.55). The term in the numerator in square 
brackets raised to the power (— 1 — e) changes sign on the integration path. It means that 
the corresponding integral has an imaginary contribution. We can rewrite 11x^x2 follows: 

II^,l2 = ^^-^ -f . (4-60) 

(1 + ^)5(1-0:1(1-1(1-0:2))) 

The integration of 11x^x2 ^^^^ ^1 more difficult because of the additional term Xi^ . We 
proceed by expanding Xi^ as (1 — elnxi + yln^xi + ...). One can see that only the ffist 
term of this expansion gives rise to a divergence in the subsequent integration. As we will 
need to find a subtraction term for this term we will treat it separately. For the remaining 
terms we do an overall e-expansion of (4.60) and then perform the remaining integrations. 
We therefore substitute x'^^ by 1. The integral looks simpler (we denote it by 11^^^^). 
Integration over xi yields 

_ {sx2)-^-^X2{-l - ^S)-^ 2F^il, -e, 1 - e, -j^S^W^)^ 



^^"^ e(l+£)(l-X2)(l + (f:-l)x2) 

[SX2) X2[l~X2) 2^l(l,-g, 1 -g:- ) 

s{l + e){l - X2){1 + {i- 1)X2){§X2)-' 



(4.61) 



Note that we omit the imaginary shifts i5 in the arguments of the hypcrgcometric functions 
2F1, as the branch cuts of 2-F1 are never crossed in the physical region. If we would directly 
integrate the above expression we would have a divergence at ,X2 1. We must therefore 
define a subtraction term. If one uses 11^^ in the present form the definition is rather 
difficult: as ^2 ^ 1 the argument of the first function 2-F1 goes to infinity. To circumvent 
this problem we can use one of the relations between hypergeometric functions to transform 
the argument of the function. As a result we pull out the divergent term as an overall 
factor multiplying the hypergeometric function with a transformed argument. The whole 
expression can be rewritten as 

[i--X2) i - £, z(i+(_i+i\xo) ^ 

11x2 = 7 7 X ^ ^ ' ' + (4.62) 

5 (l+(-l+t) ,T2) (l + £)5 



{§{1 - X2)r'-'xr (1 + + X2y " (-1 -iS)-^2F^{-e, -e, 1 - e, ^^r^y^jC ) 

(^(1 - X2f + s X2y\l + e) e 

Now all the poles arise from the factors (1 — 2:2)"^"^^ and (1 — 2:2)""'^"^, and we can derive 
the necessary subtraction term following the above procedure. We briefiy mention that 
when X2 — 1 the first hypergeometric function (in the first line of Eq. (4.62)) takes the 
value 2F1 — 1. The second hypergeometric function takes the value —en/ sin(— stt) which, 
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in turn, can be expanded to order £^ as l + e'^({2) + 7£^C(4)/4. Thus, the subtraction term 
reads 



,,0. (1 - ^2)-^^^- , (1 + ^^C(2) + ^) (1 - x,)-^- (-1 - ^6)-^ 

{l + e)esi {l + e)ss ■ K ■ ) 

Integrating this subtraction term we arrive at 

which can finally be expanded up to . 

Now that we have found a suitable subtraction term we can proceed with the remaining 
terms. We subtract from 7/°^ (Eq. (4.62)) the subtraction term 11^'^ (Eq. (4.63)). Since 
the result is convergent with regard to the integration over we can expand the result 
in terms of e before integration which greatly simplifies the problem. Wc then do the last 
integration. The difference //^^ — 11^'^ must be expanded up to third order in e. The 
reason for this is that we have already one pole ~ 1/e after the Xi-integration. Therefore, 
in order to get results up to second order the hyper geometric functions have to be expanded 
to third order. The expansion for one of the hypergeometric functions is done using (4.57). 
For the £-expansion of the second hypergeometric function one gets 

2Fi(-e, -e, 1 - e, ^) = 1 + \^i<,{z) (4.65) 
-e^ ln2(l - ^) ln;s + ln(l - ^)Li2(l - z) - - z) - Uj,{z) + C(3)' 

Using these results for the s-expansions we expand //^^ — //°2* up to and integrate the 
resulting expression. Finally, carefully collecting all the relevant pieces, we arrive at the 
final result for our second four-point function. In order to reduce the length of the final 
result for D2 we introduce four more abbreviations. We write 

D = m^s — tu^ /d = In — --1 = In — ^, Z[/ = In — - (4.66) 

mr 

Our result for the four-point box diagram D2 reads: 
ReD^-') = 2/(st), ReD^-') = -[/,+2y/(st), (4.67) 



ReL>f = [2Z,/,-5C(2)]/(st), 



ReD 



(1) 



1 

71 



^ + 3/? + ^ - ll (4/t + 3/. + 4/,3 - /.4) + \ l%4 - IJl, + I 4- 



^ll {6lt + 3/, + 4/,3 + 2U) - Is - \ii - Uz^i - 2lt + 2/,3 + /m)) - 
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m 



Ix - C(2) + C(3) + 2 Li2 — + Li2 



^5 



2m2 



(^.-2Zt)- 



4Li2 — T Mt - 2Li 



- 2Li2 - Ut - 4Li3 — + 2hh{-x) - 



2Li3 

2Lis 
1 

s7 



2Li, 



2^5 



+ 2Li3 



/ m 

2^6 



^m2(l + /3)^ 



4Li3 ( 4 



-2Li3 



2^5 



2^6 



23 2 2''' 3 '''2 2 212 

-^Ig^t^x + ^hh^x ~ -^I't^x ~ -^s^tIx ~ ^Ishlrlx ~ 3/j/r^a; + 3/t/y/a; — -IJu^x + 

A 1 1 l2, ^ 7 7 72 72 72 ^ 7 7 72 ^ 7 7 72 7 7 72 ^ 7 73 7 73 



12^^ 16" 24 



19 



2 2 2 5 2 -'-1 

ltlxlz3 — ^^sh^xhs ~ iOk^T^xhs + -^h^lhs + '^k^lhs ~ ^^T^x^zS ~ 24^^x^^3 ~ 



5 I XI 



1 

2' 



I' 



2I22I3 3 3233 2 2 

16 5 

IsItIuIz4 + hlu^zi — hlu^zA — Irluh^ + '^tyxlzi ~ ^Ishlxh^ + ^ll^x^zA + ^Idrlxlzi 

3 5 1 

—IglulxlzA: + hlulxlzA — "Hrlulxlzi + -h^x^zA + IJx^zA + -It^x^zA + -^K^zA + 

Z Z 6 

^IglzshA — ^^ilglflzslzA — ^h^z^^A + ^h^rh'ihA + ^U^rh'ihA — ^IglJ'z'ilzA + 
^kluhslzA — 2lslxlz3lzA — ^ItlxhslzA + ^IxlxhslzA — ^il^lzslzA + ^h^l^lzA — 
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2 22 2^22 2 2 2 2 

15 2 1 

3 1 

'^Ishlz^.lfS + 2/s/x-^z3^/3 + '^klxhs^d + 2/344/a + 2lslxlzdp — '^hhshd/S ~ 
'^Ixhshilp ~ 2/^/^ — /s^t^^ — ^Iglxl'^ — hl'xl''^ + I'slzsl'^ + Ixhz^'^p + hh^lp + 

5ll — 61tIz3 ~ '^Ixhs ~ 6Z^3 ~ 2ZrZ^4 + 2Zj;Z;j4 — 12lz^lzi — 71^4 + hi^lr + 21Zx + 

16/^3 + 6/^) + 4M/3 - 4Za:i/3 - 4:lz3h + hi-'i-'^k - 4/t - + 2/^3 + 10/^4 + 

2/^) j C(2) + (7/, - 8lt + 4/r - 3/, + 2/,3 + 4/,4)C(3) - yC(4) + 4Li2 
Zt^x - 2Li2(a;) (l^ - Ut + Uzs - Ut + '^hh - Uzs) + 2Li^(-x) + 
Li2 + 2ls{lt - Ix) + '2lxlt + it) + Li2 (-211 - ^kh - 2Ut+ 

ll + 2k {2lt -It- Ix)) + Li2 (^) (-2Li2 (-^) - 2Li2(-a;) - y/^ - 2lTk 
-2ll - 4 - 2lzslz4 - 4 - kh + eitihs + - 14C(2)) + Li2 j {^ll- 
ll - 2hi3lt -It + In) + 2/t(3/r + Q + 2(/t - /x)/. - 8C(2)) + 
Li2 (--) (3^,2 - 2k{?,lt - 21t + lu)-2 (ij + 2ll - kiAh + k) + 1^) - 8C(2)) 



X 



+Li2 ( -4t) f-2Li2 ) - 4/3/* - 2ll + - - 4C(2) ) - 



mH 

T 



m? 



Li2 (2Li2(-x) + ll - AIM + Mk + ^It - 2lx - lz3 - 2lz4) + 2C(2)) + 
Li2 {-ll + 2k{-k - SIt - Six + hs + 21za) - 2C(2)) + 

Li2 ( ~^2^^2^^ ) (-2Li2(-x) - 2ll - Alth - 21xIt + + 2/^(2/^ - h + lx)+ 
8C(2)) + U^i-x) _ Qij^ ^ 7^2 _ 2/^/^3 + gy^ _ 2/^/^3 + 5/^/^ _ Qij^+ 

\lslx + 2/,/,3 + 12C(2)) + Li3 [ ^ j (-6/, - + /r + 2/, - 4/,4) + 



27 



Li3 (^) (-4/. + 8lt + 3/t + 21, - 61,3 - IO/.4) + 2U3 (^^) {k + h - - 
hs - 2U) + 2 (Li3 {^'^^~/^ ) - Lis [^^^) ) i^t + SIt + SI. - - 

2U) + 2Us {k + It + 5k - lz3 - 21,4) + 2Li3 (3/, - 3lt - 21t) + 

Tn^t\ T . ( t^W , , , , N ( T ■ ( Tn?u\ ^ . f Z3 



2 ^Li3 ^-_J - Li3 j j (2(, + - iW) + 2 (^Li3 j - 2Li3 ^ _^ 

+Li3 -k-r It) - 2Li3 (- J) (/^ - /t) - 2Li3 {k - 2/r) + 

2Li3(x)(3Z, - Q + 2Li3 (/. - Z.) - 2Li3 ( ''^\~^^ ] (21, -lt-lT + 

y , ,,3 (I) , , (1)4.2 (l,3 {^^) . L.3 (I) . Li3 (I 
X (Z, + Q + 2Li3 (^3^) (-2^^ + + /t + i.) - 2Li3 j + 4/^ ■ 

5/. - 4^,3) + 2Li3 (yT^) + - 21, - 21,3) + 6Li3 (^t - - 

/,3) - 2Li3 ( 2^/1^ 1 (2^^ - ^* - 3/t + 2/, + /,3 + 2/,4) + 2Li3 ' ^ 



x(2Z, - - 3Zt + 21, + Z^3 + 2/^4) + 2Li3(-x)(8Z, + 4/^ - SIt - 7lx - 21,3 + 
4U) - - 6m (4) + 3m Q + 16m (^) - 12m 

-12Li4 + (1, 0, 0, ^) + h.^^ (1, 0, 0, -| 

L_++ (^1, 0, 0, - L_++ (^1, 0, 0, ^ j - 5L_++ (^1, 0, -1^ - 

/ rn"^ i\ ( m'^ \ 3 ^ ( rri^ rri^ \ 

2L_,, (1. 0, 3^. -] - 2L_,, (^1, 0. -„xj - (1. Zr- Zt- "] - 

( ^ u\ / l\ / ^ 

4L-.. (1, Zt' ' 7 j + 2i^-.. (1, Zt' ' j + ^i^-.. (l, z^' Z^' 

-3i^-.. (|, 0, 0, + 2L, (0, 0, 1) + 2L, (0, 0, - J, .) - 

2L+ (0, 0, y , + 2L+ (0, 1, -1, ^) + 2L+(0, 1, -1, x) + 4L+(0, 1, ^, — ^ ' 
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(4.68) 



Isik + I + /.3 + /m) - 2C(2) - 2Li2 - 2Li, (^^^^ 



~ '^^xkh^ — h^z^ + hi ~ ( 2/j — -/^ — 2/t(/a. + Iz^) — 



TT 
St 



2/:.(/.3 - Izi) - 4) - K ( 4^0. + /.3 + Iza]+ 4(/i + /:.)C(2) + 2Li2(a;) + 



2Li2 



-2m2 



^3 
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{Is - Ix) + 2Li2 - + 2Li2 — + - 2Li3(x) 



^5 




^2 ^4 

FIG. 4. Massive box with two massive propagators. 

The diagram corresponding to the third four-point function with two massive prop- 
agators is shown in Fig. 4 where we write 



D2, = D{-p2,p4, -pi,0,0,m,m) . 
The kernel (4.40) for can be written as 

K = act + bdu + (c + d^, 
where we have introduced the positive- valued dimensionless variable 

u> 1, s > u. 



u 



u 



(4.69) 
(4.70) 
(4.71) 



For the Feynman parameterization we choose {a,b,c,d} as {xi{l — X2), 1 — Xi, a;ia;2(l — 
Xs), X1X2X3}, which gives the following integrand: 



xlx2 



XiX2 + tXiX2{l — X2){1 — X3) + UXi{l — Xi)x2Xs 



-2-e 
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The above expression never becomes negative. Therefore, the entire result for the box D3 
does not have an imaginary part. One can set 5 = in the kernel from the very beginning. 

That the box posesses no imaginary part can be seen in a less technical way by 
appealing to the Landau-Cutkosky cutting rules. The diagram corresponding to the box 
D3 shown in Fig. 4 does not admit of any cuts such that the cut hues of the diagram are 
on their mass shell simultaneously. 

As before we obtain two terms after the first integration over x^. They are 



-l-2e -l-e 



X2+i{l-X2j) 



-1-6 



(1 + e)[u{l - xi) - txi{l - X2)] 

jjDs ^ X{'X2^'' [XiX2 + m(1 - 

"^"^ ~ (1 + e)[u{l - xi) - ixi{l - X2)] ■ 



(4.72) 
(4.73) 



Note that the denominators of I^^^^^ and IIxiX2 change sign on the integration path, while 
the numerators stay positive (i.e. the relevant integrals have branch cuts). This can easily 
be seen by considering the numerators and denominators of the above integrands at two 
particular values of the variable xi, for instance at xi = and xi — 1. This means that 
although the whole box integral does not have an imaginary part, the two terms in (4.72) 
and (4.73) separately give rise to unphysical, spurious imaginary contributions. Of course, 
these are artefacts of having split the result into two terms. On the one hand, this somehow 
complicates things. On the other hand, the cancellation of imaginary contributions in the 
sum of the two terms (4.72) and (4.73) will serve as a good check for our final result. To 
control the imaginary contributions we do the following replacement in the denominators 
in (4.72) and (4.73): 

u{l — xi) — ixi{l — X2) — > u{l — xi) — TXi{l — X2), T = i — iS. 

We start with the a;2-integration of the term /^^^ Eq. (4.72): 



X2 (x2 + X2))] ^ ' 2i^i(l, -2£, 1 - 2£, 



U+T—X2T 



) 



2u{l + e)e 



(4.74) 



The above expression is singular at the lower integration limit X2 — due to the term X2^~^- 
To find a subtraction term, we follow exactly the procedures defined after Eq. (4.56). Our 
subtraction term reads 



jD3,S _ 



X2^-' t-'-' [-1 - 2e In (-1) + 4£2Li2 (^) + 8e' Lig (^) + We^Lu 



2u{l + e)e 

The above subtraction term can easily be integrated over X2 to obtain 

r^-^ [1 + 2e In (-1) - Ae^U2 {^^) - 8e' Lis (^) - l^e^U, 



2{l + e)e^u 



(4.75) 



(4.76) 



31 



which can readily be expanded in £. As the difference I^^ — 7^^'* does not contain any 
poles, we can expand the difference in a series in e and perform the analytical integration 
over the last variable X2. 

To integrate the term (4.73), we split H^ix^ contributions, i.e. 



X2 ^ [xiX2 + u{l — Xi)\ 



-l-£ 



X1X2 



(1 + s)[u(l - xi) - rxi(l - X2)] 

(xr^ - 1)X2^~' [xiX2 + U{1 - Xi)]~^~' 

{1 + e)[u{l - xi) - rxi{l - X2)] 



(4.77) 



Then we integrate the first term in (4.77) over xi to obtain two hypergeometric functions 
which are expanded up to As was done previously, we then introduce a subtraction 
term similar to (4.75) which is integrated analytically. The finite difference of the original 
integral and the subtraction term is then ready to be integrated over the last integration 
variable. 

For the second term in (4.77) we introduce a subtraction term before the last two 
integrations: 



D3,s 



{x^'-i)x2'-' m-xi)] 



-1-e 



{1 + e)[u{l - xi) - TXi] 
which is obtained from the second term in (4.77) by the substitution X2 



(4.78) 



in all terms 



except for x. 



-l-e 



We first trivially integrate out X2 in (4.78) and expand the resulting 



expression in a series of e. Because of the factor (x^^ — 1) this expansion starts at the order 
£°. Thus, the subtraction term is finite and ready for the last integration. 

As the difference of the second term in (4.77) and its subtraction term (4.78) does not 
contain any poles, we expand it up to e'^ and integrate over Xi. Finally, collecting all the 
relevant pieces, we perform the last integration. Our result for the four-point box diagram 
D3 reads: 



ReDf = 



l/{tu), 



-[lt + lu]/{tu), 



(4.79) 
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The non-planar topological structure of the four-point function D2, implies that has to 
be {t <-> M)-symmetric (see Fig. 4). This can best be seen by exchanging the momenta 
P3 ^ P4 in Fig. 4 followed by a twist of the r.h.s. of Fig. 4 ^. The {t <-> ii)-symmctry 
provides for a check on our results for D3. The symmetry is obviously satisfied for Y{.eD^^\ 



ReDs"^"* and ReDg^'' but is not manifest for YieD^C and ReDg^'' in (4.79). However, it is 
quite straightforward to verify numerically that the (t<-^ti) -symmetry indeed holds for all 
coefficient functions in (4.79). 

Apart from the internal checks mentioned earlier the most important check on our 

^The {t u)-symmetry is not so easy to see when exchanging pi ^ p2 in Fig. 4. In this case the 
(t<-^u)-symmetry becomes apparent only after Feynman parametrization. 
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four-point function results has been a comparison with numerical results provided to us 
by M.M. Weber [17] for several phase space points. Within numerical errors we have 
found complete agreement with the results of M.M. Weber for each of the three four-point 
functions. It is important to emphasize that the approach of M.M. Weber to numerically 
evaluate the four-point functions is completely different from ours [18]. 

V. SUMMARY AND CONCLUSIONS 

We have presented analytical results up to 0{e'^) for all the massive scalar one-loop 
integrals that arise in the calculation of one-loop matrix elements in heavy flavor hadropro- 
duction. Many of our results are new (sec Table 1). The one-loop scalar integrals are 
needed for that part of the NNLO hadroproduction of heavy flavours which is obtained 
from the product of one-loop contributions called loop-by-loop contribution. 

What remains to be done in order to obtain the full one-loop amplitude structure is to 
take into account positive powers of e (up to 0{e'^)) resulting from the Passarino-Veltman 
decomposition and the spin algebra. The full one-loop amplitudes to order were given 
in [13]. The missing results for the e- and e^-coefflcients of the one-loop amplitudes will be 
presented in a forthcoming publication [19]. In a last step the amplitudes themselves have 
to be squared, which, in the case of gluon-initiated production, will generate further positive 
powers of e in dimensional regularization. The calculation of the loop-by-loop contributions 
in Fig. 1 is a necessary starting point in the evaluation of the NNLO contributions to heavy 
quark pair production in hadronic interactions. It is very likely that the calculation of the 
other three classes of diagrams in Fig. 1 will proof to be very difficult. This holds true in 
particular for the massive two-loop box contributions. 

In the Laurent series expansion of the scalar one-loop integrals the successive coefficient 
functions increase in length and complexity with each order of e. The reason is that the s- 
expansion of the integrand before the last parametric integration itself generates coefficient 
functions with increasing complexity with each order of e. The most complex expressions 
arise from the box contributions where one encounters multiple polylogarithms up to weight 
and depth four at 0{e'^). 

In a numerical NNLO evaluation of heavy hadron production the various contributing 
pieces will have to be evaluated at many values of the kinematical variables. This requires 
efficiency in the numerical codes for each of the contributing pieces. We believe that we have 
provided for such numerical efficiency in the loop-by-loop portion of the NNLO calculation 
by presenting results in analytical form which are fast to evaluate numerically. All our 
results are available in electronic form [20]. We are planning to present our results in terms 
of multiple polylogarithms in the near future. In recent years number of new methods were 
developed for semi-numerical evaluation of general Feynman diagrams (see e.g. [18,21,22]). 
First numerical tests have shown that our representation in terms of the L- functions perform 
better than the present implementation of the flexible all-purpose approach described in 
[18,21]. 
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The analytical results presented in this paper cover the whole kinematical domain with 
a single expression. They evaluate numerically very fast and efficiently. Further advantages 
of having the results in analytical form are that they allow one to investigate various limiting 
cases as well as their analyticity properties. Also, when analytical results are available the 
mathematical structure of the results becomes manifest which would not be visible in a 
purely numerical approach. 

The full calculation of the NNLO corrections to heavy hadron production at hadron 
colliders will be a very difficult task to complete. It involves the calculation of very many 
Feynman diagrams of many different topologies. The problem is further complicated by 
the fact that heavy hadron production is a multi-scale problem with three mass scales 
provided by the kinematic variables s and t in the loop expressions, and the mass of the 
heavy quark. It is clear that an undertaking of this dimension will have to involve many 
theorists and cannot be done by a single group alone. In this sense the present calculation 
is a first step (or second step [23]) in the direction of obtaining NNLO results on heavy 
hadron production at hadron coUiders. The present calculation allows one to obtain a first 
glimpse of the mathematical and computational complexity that is waiting for us in the full 
NNLO calculation. This complexity does in fact ah^cady reveal itself in terms of a very rich 
polylogarithmic and multiple polylogarithmic structure of the Laurent series expansion of 
the scalar one-loop integrals as shown in this paper. 
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In this Appendix we evaluate a special two-point integral which is needed for the cal- 
culation of the one-loop fermion self-energy diagram insertion into the massive external 
fermion line. This integral is also needed for the definitions of the fermion mass and wave 
function renormalization constants in the on-shell renormalization scheme. In particular, 
we need to evaluate the integral 
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up to 0{p'^ — m^). We therefore Taylor expand h around — w?: 



h 



h 



+ 



dh 



= Eq + Ei{p'^ - m'^) + . . . 



{p^ — m^) + . . . 



(A2) 



Note that the expansion coefficients Ei in (A2) are functions of e. The ffist coefficient Eq 
is nothing but obtained in Section II. The second coefficient Ei is proportional to the 
sum of a scalar and a vector integral obtained by differentiating 7i w.r.t. p^. One obtains 



El 



where 



One finally has 
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e{l - 2e) 
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+ 0[{p'-m'y]. 



(A3) 



(A4) 



(A5) 



The result for the integral Ji in the form (A5) was used in [24] to evaluate external heavy 
quark self-energy diagrams and obtain heavy quark wave function renormalization constants 
in the NLO calculation. 



APPENDIX B 



In this Appendix we shall demonstrate how the L-functions introduced in Eqs. (3.20) 
and (3.21) are related to multiple polylogarithms as defined in [14]. Multiple polylogarithms 
are defined as a limit of Z-sums, e.g. 



_1 -^2 ■■■•^k 



The number w = mi + ... + mfc is called the weight and k is called the depth of the multiple 
polylogarithm. The power series (Bl) is convergent for \xi\ < 1, and can be analytically 
continued via the iterated integral representation: 

xiX2...xi^ / 7, \ mi— 1 ,, 

dt \ dt 



'dt \"''~' dt (dt \ dt 



-o o ... o — o 



t I X3...Xk — t \t J 1 — 



X2X3...Xk - t 

(B2) 
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where the following notation is used for the iterated integrals: 

I ^o...oJ^= j I ^^"-^ x...x/^. (B3) 

an — t ai — t dn — tn O^n-l ~ tn-1 ^ Oi — ti 

Note that the classical polylogarithms 

Un{z) = J ^^^^^j^dC, n>2; Lh{z) = - ln(l - z) (B4) 



are a subset of multiple polylogarithms. Examples of this statement can be found in the 
subsequent discussion. 

We start by considering the single-index L-function Eq. (3.21): 



lii(«i + (Jiy) 

Lai{ai,a2,a3,a4) = dy ■ Li2(q;2 + as?/)- (B5) 

Jo aA+y 



™i ln (ai + (Ti y)^._ 

After changing the integration variable y — {ti — 0:2) / as we get 

'''T'dU Infai + ai'-^) "^1"' In ^ + ln(™ - 0^2 + ^i) 

^ CKs q;4 + — — - J a^a^ — q;2 + ii 

OC2 '^3 a2 

The integration interval can be split into two pieces, [q;2, 0] and [0, a2 + "a]. One can then 
write as a sum of four terms: 

where we have introduced the notation 

a = 7 = a^a^ — a2- (Bsj 

Regarding Eq. (B7) it is clear that there are only two different types of integrals to be dealt 
with: 

/^Li2(t0 and y dU^^^^^U2{U). (B9) 

with upper limits = 0:2 + "s or tm = (^2- The first integral can be evaluated analytically 
in terms of standard logarithms and classical polylogarithms up to Lis. However, the 
same integral can also be expressed in terms of multiple polylogarithms via the integral 
representation (B2), e.g. 

^7 + ^1 J 'J + ti J t2 J I - ts V 7/ 

^ 
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We now deal with the second integral in (B9). Consider the following multiple polyloga- 
rithm of weight four: 

a tm\ f dt2 f dti -r • / \ f dti f dt2 



Li2,i,i -7, — = / r / r^i2{ti) = / Li2 ti) / — 

V -aj J -a-t2J -^-ti ^ + *i J a + t2 



/(It f dt 

—^U2{h)H(^ + tm)- / —^U2{h)hi{a + ti). (Bll) 



In the first step we have used the usual trick to change the order of integration. As already 
noted before (see Eq. (BIO)) the first term on the second line can be expressed through a 
multiple polylogarithm of weight three. Thus one has 

Jdt,^^^^^U2{t,) = -Lz2,i,i (-7, ^, ^) - L^2,^ (-7, ^) ln(a + t^). (B12) 

Finally, substituting Eqs. (BIO) and (B12) into Eq. (B7) we arrive at the desired relation 

L^^{ai, a2, 0:3, Qi4) = -^'^2,1,1 I -J, ^, 3^ j - -^'^2,1,1 f "7, ^, -^Z^ 1 ^^^^^ 



+Li2A i -7, — ^ I ( In — + ln(Q; + a2) 
V 7 / V as 

-Li2A I -7, — I ( In — + ln(Q; + a2 + 0:3) 

\ 7 / V CKq 



r f \ ^ ln(Q;i + (7iy) ln(Q;2 + (72y) ln(Q;3 + (73y) , . 

Va2<T3(«i,"2,tt3,tt4) = / dy ■ (B14) 

Jo OiA+y 



1 J \ 012, 
7 y V q;3 

where a and 7 are defined in Eq. (B8). 

Next we turn to the triple-index L-function Eq. (3.20) 

"i ^^}n{ai + aiy) \n {a2 + a2 y) ln(Q;3 + asy) 
a4 + y 

We again change the integration variable y — {1 — — ti) / and obtain 

- - ln(7i - ^h) ln(72 - ^t,) Ml^, (B15) 

J as (73 CTs 73 - ti/as 

1— as 

where 

71 = + — (1 - «3), 72 = 0:2 + — (1 - CKs), 73 = 0/4+^^ —. (B16) 

(T3 (73 (73 

At this point we take ai = (72 = —(73 = 1, which does not affect the generality of our 
further discussion. Splitting the integration interval in (B15) into two pieces as before, we 
arrive at 



1— as— (T3 1— as " 



L++_ = ( / - J ]dti ln(7; + h) ln(7^ + t,) ^^'l (B17) 



73 + ^1 
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with 

7i = «! - 1 + Q;3, 72 = «2 - 1 + as, Ts = «4 - 1 + ^s- (B18) 
Now consider the following multiple poly logarithm of weight four: 




where we have introduced the notation 




(B20) 



The third term on the last line of (B19) is the integral of the required type needed to 
express !/++_ in Eq. (B17). We can write: 

7dt, ln(7i+ti) ln(72+ti)i^^^^^ = Li^,,,^,, f-73, ^, ^, ^] -I\tJ+I"{tJ. (B21) 
^ 73 + 11 V 73 72 71 / 

To express the last two terms of (B21) in terms of multiple polylogarithms we proceed as 
follows. Consider first 




(B22) 



40 



from which we immediately conclude that 

/"(i^)/ln(7i + i^) = Lzi,i,i f-73, ^, ^] + In (7^ + tj Li,,, f-73, ^\ . (B23) 

V 73 72 / V 73 / 

Note that the above multiple polylogarithms of weight three and weight two can also be 
expressed in terms of logarithms and classical polylogarithms by direct evaluation of the 
corresponding integrals. 

For the first integral in (B20) we write 

/'«„) = = Ut}<!l±}Au,, U, ^] . (B24) 

{ ii + u I 73 + ^2 I ii + u V 73 y 

On the other hand one has 



Li.,,, f-73. 2i. ^, ^) = / ^ / -^U,, (-,3. ^) = (B25) 
V 73 7i 72 y ^ 72 + ^2 ^ 7i + ^1 V 73 y 

^ 7i + V 73 y I 72 + 12 
ln(72 + t^) / —rLii,i -73, - / dti — — Li^^i -73, = 



-ln(72 + t„i)Lii,i,i (-73,—,^^) -I'{tm)- 
\ 1^ 1\ J 



We then conclude that 



I\tm) = -Li,,,,,,, ( -73, -, -, - ln(72 + tr,,)U,,,,, ( -73, ^, ^ ) . (B26) 

V 73 7i 72 y V 73 7i y 

Finally, substituting Eqs. (B23) and (B26) into Eq. (B21) we write down the result for the 
integral 

/ dt, ln(7i + t,) ln(72 + ^ ^^^^^ ^ ^^^^ ^ tjLi,,, f-73, ^] + (B27) 

^ 73 + 11 V 73 y 

ln(72 + tm)Li,,i,, ( -73, li^ I + in('y^ + im)-^«i,i,i ( -73, — , ) + 
V 73 7i y V 73 72 y 

r- / 72 71 . / 7l 72 -^m^ 

Lt,,,,,,, -73, — , — , + Li,,,,,,, -73, — , — , 



73 72 7i y V 73 7i 72 y 

which demonstrates how the two integrals in Eq. (B17) representing L^^_ are related to 
multiple polylogarithms. The discussion of the other cases for the Lfj^^^^^ functions proceeds 
along similar lines. 
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APPENDIX C 



In Appendix C we discuss properties and identities involving the single- and triple- 
index L-functions in Eqs. (3.20) and (3.21). There are two different categories of identities 
which we discuss in turn. We start by considering the simplest identities originating from 
symmetries related to permutations in the indices and arguments. We then present further 
identities based on integration-by-parts techniques. 

1. Symmetry properties 

We start with the single-index function L^j^{ai,a2,as,a4). One notices that a change 
of the integration variable y — > 1 — y results in the identity 

Lo-i(q;i, q;2, as, 04) = --L_<^-^(ai + cti, 02 + as, -Qia, -"4 - 1) (CI) 

which implies that L_ can always be related to L+, and vice versa. We have thus written 
our results for the three-point and four-point functions in the main text only in terms of 
the L+ functions. 

Next we turn to the triple- index L-function. Note that Lo.^o-2o-3(q;i, Q;2, eta, q;4) is sym- 
metric under permutations of any two pairs of indices and arguments {uj, cti} and {aj, aj} 
for (i 7^ j). The same change of variables as above y — > 1 — y results in 

Laia2<T3{0il, Oil-, as, ^4) = --L_o-i-(T2-o-3(ai + ^2 + 02-, ^S + (^3, -OLi, - 1). (C2) 

Therefore, from the eight functions L , L |_, !/_+_, L_| , !/_++, !/+_+, !/++_, 

only two are independent. We have chosen to write our results in terms of and 1^+++. 

2. Integration-by-p£irts identities 

The triple- and single-index L-functions L_|__|__)_, L__|__|_ and defined in Eqs. (3.20) 
and (3.21) have been devised such that they have neither branch cuts nor poles on the 
integration path y G [0,1]. This also implies that the and L_|_ functions are 

real. Remember that the branch cuts for the In and Li2 functions are (— oo, 0] and (1, +oo), 
respectively. The domains of the functions L_+_|_ and L+ are 

L+++(q;i, 0:2, cts, ^4) : tti > 0,0:2 > 0,0:3 > 0, 04 < —1 or 04 > 0; 
L_++(oi, 02, 03, 04) : oi > 1, 02 > 0, 03 > 0, 04 < — 1 or 04 > 0; (C3) 
L_|_(oi, 02, 03, 04) : Oi > 0, 02 < 1, 02 + 03 < 1, 03 7^ 0, 04 < — 1 or 04 > 0. 

Looking at the definition of the triple-index L-function in (3.20) one concludes that the 
boundary points Oi = and/or 02 = and/or 03 = can be included in the domain of 
definition for L_|__|__|_. The same holds true for oi = 1 and/or 02 = and/or os = for 
L_++. Also, from the definition of the single-index function L+ in (3.21) on concludes that 
the boundary point ai = can be added to its domain of definition. 
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The points q;4 = {—1,0} can also be included in the domain if the values taken by 
other parameters aj guarantee the convergence of the integral. In what follows we assume 
everywhere in this appendix that the conditions (C3) are satisfied. Nevertheless, it is always 
possible to analytically continue the parameters to the complex plane. 

In order to obtain integration-by-parts identities one makes use of the standard integration- 
by-parts formula 



J UV'dy = UV\^^ - J VU'dy . 



(C4) 



We start with the triple-index functions and defined in Eq. (3.20). Set- 

ting U equal to the numerator [ln(Q;i -|- aiy) ln{a2 + (T2y) ln(Q;3 -|- a^y)] and V equal to the 
remainder {a^ + y)~^ we then arrive at 



-L+++(Q;i,Q;2,Q;3,Qi4) 



Q!4 > 



q;4 < —1 



ln(ai + y) \n{a2 + y) ln(a3 + y) ln(a4 + y) 
-L+++(q;4, q;2, «3, cti) - L+++(q;i, a^, as, ^2) 

ln(ai + y) ln(a2 + y) ln(a3 + y) ln(-a4 - v) ^ 
-L_++(-q;4, 0:2, a;3, ai) - -L_++(-q;4, ai, 0:3, 0:2) 
-L_++(-q;4, ai, a2, as); 



(C5) 



and 



-++ 



[ai, a2, as, "4) 



a4 > 



a4 < —1 



(C6) 



ln(ai - y) ln(a2 + y) ln(as + y) ln(a4 + y) ^ 

-L+++(a4, a2, as, -ai) - -^-++(a;i, Q;4, ^3, ^2) 
-L_++(ai,a2,a4,a3); 

ln(ai - y) ln(a2 + ?/) ln(as + y) ln(-a4 - 2/) ^ 
-L_++(-a4, a2, as, -ai) 
+i^-++(as -t- 1, ai — 1, — a4 — 1, — a2 — 1) 
+I/-++(a2 -I- 1, ai — 1, — a4 — 1, —as — 1). 

For the second part of Eq.(C6) we have made use of relation (C2). 

There are some special cases when some of the a^ take values on the boundary of the 
domain of definition where one can still make use of the identities (C5) and (C6) even if the 
conditions (C3) are not met. For example, for the case {ai = 0, a4 = —1} the identitity 
(C6) is still valid. There are similar special cases for further identities to be derived below. 

The intcgration-by-parts identities for the single-index function are more involved. 
To prepare ourselves we first write down the derivative of the dilog function in the integrand 
of (3.21). One has 

dLi2(a2 + a^y) _ ln(l - a2 - asy) 



dy 



— + y 



(C7) 
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In the case of the single-index function it will prove important to consider two different 

choices for U. Wc start by setting the whole numerator [ln(Q;i + aiy)Li2{a2 + Oi^y)] in the 
integrand of Eq. (3.21) to U. For V we then have (0:4 + y)~^. One obtains 



a4 > 0,as > 



q;4 > 0, aa < 



L+(q;i, a2, as, 0:4) = 

ln(Q;i + y)U2{a2 + a^y) ln(a4 + v) ^ 
-L+(q;4, q;2, as, a^) 

, ai, a4, + ln(«3) / '^^^^^^^^^^dy- 

1 




q;4 < —1, 0:3 > 



ln(ai + |/)Li2(a2 + a^y) h\{ai + y) 
-L+(q;4, q;2, "3, tti) 

+^+++(^' «4, + ln(-«3) / '"("^+gM"^+^) ciy; 

1 





ln(Q;i + y)Li2(Q;2 + as?/) ln(-Q;4 - y) 
+L^{—ai — 1, 02 + eta, —eta, — ai — 
-L_++(ai + 1, -a4 - 1, - 



(C8) 



0:2 +^3 

as 



) 



^4 < — 1, q;3 < 



ln(ai + y)lj\2{a2 + a^y) ln(-a4 - y) 
+L+(— 0:4 — 1, 0:2 + 0:3, —0:3, —ai — 

+L_++(-a4, ai, + ln(-a3) / '^^^^^^^^^dy. 



An additional condition for (C8) has to be explicated because it does not follow automati- 
cally from (C3), namely the parameters a2 and are restricted by 



— < -1 or — > 0. 



aa 



aa 



(C9) 



The integrals in (C8) are simple enough to be evaluated in terms of classical polylogarithms 
up to Li3. We do not provide explicit results for these integrations since they are rather 
lengthy and, in addition, depend on relations between the parameters. 

A second choice for U in (C4) provides further identities for L^. In this case one sets 
112(0^2 + O-zy) to U and ln(Q;i + y)/{(y.A + y) to V . To calculate V one has to differentiate 
between three cases for the set of parameters ai and q;4. One has 



V 



I 



\n{ai y) 
a4 + y 



dy 



Q.l < «4 
ttl > 0:4 > 
^4 < — 1 



Li2( 



ln(«i + y) ln(^: 
ln(ai-«4)ln(a4 + y)-Li2(^); 
ln(Q;i — a^) ln(— 0:4 — y) — Li2( 



(CIO) 



Using Eq. (CIO) one then obtains 



L+(q;i, q;2, q;3, 0.4) 
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cti < 0:4, 0:3 > 



U,{a, + a,y) {\n{a, + y) In(^) + U,{^^ 



+L_++(i^,ai,«4 ^ 



1 ln(Q3) 



1—02—03 



03 ' 



02+03. 



03 



«3 



) + 



-ln(a4-ai) ln{ai+y) ln{^-^-y) 



dy; 



ai < a4, 0.3 < : ^2(0^2 + a^y) {Ha, + y) In(^) + 



+L+++(^,ai,a4 ^ 



' 02-1 



01 



03 ' 
1 02 ' 



1 ln(-a3) 




' ai—a^^ oi— 04' 03 - 



Li2( 



Oil ~<^4 



)+ln(ai+j/) ln( 



+ 

a4 + y ■ 
CK4 — a]^ - 



^2 I 
— +V 

<^3 



-ln(a4-ai) ln{ai+y) ln(^^- — \-y) 



CKi > Q!4 > 0, Q!3 > : Li2(Q;2 + Oi^y) (ln{ai — \n{a4 + y) — Li2 

/"i- 02— 03 04+1 1 02+03 ^ I 



+ 



03 ' 04— oi ' 01—04 ' 03 
1 -ln(o3)Li2( "-^^J )+ln(oi-04) ln(o4+y) ln(l-02-03y) 

J mnf^ dy; 







a4<-l,a3>0 : Li2(Q;2 + ctsy) (ln(Q;i - q;4) ln(-ci;4 - y) - Li2( 



04 +^; 



iT / 1—02—03 04+1 1 
+ V 013 ' 04— 01 ' 01— 04' 



02+03 ' 



04— 01 



+ / 





03 ' 04— 01 ' 01— 04 ' 03 ' 

1 -ln(o3)Li2( )+ln(ai-04) ln(-04-J/) ln(l-02-a3y) 



:+y 



dy; 



ai > ^4 > 0,03 < : Li2(Q;2 + tt3y) (ln(ai -a4)ln(a4 + Z/) -Li2(^^ 



03 ' 04—01 ' 04—01 ' 03 ' 



1 -ln(-03)Li2(^l^)+ln(oi-04)ln(o4+2/) ln(l-02 -032/) 



:+2/ 



a4<-l,a3<0 : ^2(0:2 + ccsy) (ln(Q;i - q;4) ln(-Q;4 - y) - Li2( 



-LA 



a2— 1 



— —) + 

— 01 ' OT ^ 



04 +^; 
04— 01 



03 ' 04—01 ' 04—01 ' 03 ' 
1 -ln(-a3)Li2(^;7^)+ln(oi-04)ln(-a4-J/)ln(l-02-03J/) 

J dy. 

"3^'' 



(Cll) 



In deriving (Cll) it is important to take into account condition (C9). As was the case in 
Eq. (C8) the integrals in (Cll) can be evaluated in terms of classical polylogarithms up to 
Lis- 

There is also one special case of the last identity (Cll) when the first and fourth 
arguments of the single-index L+ function are equal, e.g. a, — a^. In this case can be 
expressed only in terms of the functions or as follows: 



L+(q;i,q;2,q;3,Q!i) 
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as > 



0:3 < 



\ ln2(«i + y)U2{a2 + a^y) [ + iL_++(i^, a^, «i, 



\ \v?{ai + |/)Li2(Q;2 + "si/) ^ + i-^ 



(C12) 



In(-a3) r ln^(ai+y) 







The third and last identity for the L+ function is obtained from the definition (3.21) without 
making direct use of the integration-by-parts identity Eq. (C4) . Nevertheless it can still be 
called an integration-by-parts identity because it makes use of the well-known identity 



Li2(^) = C(2) - Hz) ln(l -z)- Li2(l -z), z^C 



(C13) 



which in turn is derived from the definition of the Li2-function (B4) with the help of the 
integration- by-parts identity (C4). After transforming Li2(Q;2 + OL^y) according to (CIS) 
one gets 



-L+(q;i, q;2, as, Oii) 



as > 
< q;2 < 1, q;2 + Qia < 1 



as < 
< a2 < 1, < a2 + as 



-L+(q;i, 1 - q;2, -(^2, ^4) - -L_++(^, ai, ^, 0:4) + 

/ (C(2) - In(a3)[ln(a2 + a,y) + ln(l^ - y)]) dy; 



(C14) 



-L+{ai, 1 - a2, -a2, ^4) - -L-++(-^, ai, a4)+ 



/ (C(2) - In(-a3)[ln(a2 + a^y) + ln(^ + y)]) dy. 



A few final comments are appropiate. In spite of the rather complicated appearance of the 
identities (C5), (C6), (C8), (Cll), (C14) these turn out to be very useful to reduce the 
length of the results presented in the main text. The first step in the chain of reductions 
is to write everything in terms of the functions L+, and In a second step 

one uses the identities written down in this Appendix to find the set of arguments of L- 
functions for which the number of the functions L^, and is minimal. We 

have devised several programs for the MATHEMATICA computer algebra system which 
automatically find minimal sets of the single- and triple-index L-functions. With the help 
of these programs we have been able to greatly reduce the number of L-functions appearing 
in our results and have thereby greatly reduced their length. The same was done for the 
logarithms and classical polylogarithms using standard one- and two- variable identities for 
classical polylogarithms as given e.g. in [25]. 
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